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ABSTRACT: We investigate stability of the Higgs effective potential in curved spacetime. 
To this end, we consider the gauge-less top-Higgs sector with an additional scalar field. 
Explicit form of the terms proportional to the squares of the Ricci scalar, the Ricci tensor 
and the Riemann tensor that arise at the one-loop level in the effective action has been 
determined. We have investigated the influence of these terms on the stability of the scalar 
effective potential. The result depends on background geometry. In general, the potential 
becomes modified both in the region of the electroweak minimum and in the region of large 
held strength. 
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1 Introduction 

The issue of the stability of the Higgs potential in a flat spacetime (often under the assump¬ 
tion of no new physics up to the Planck scale) has been considered in many papers, see for 
example [1-7] and the references therein. However, the instability may affect cosmological 
evolution of the Universe and to take it into account one should couple the Standard Model 
(SM) Lagrangian to gravitational background. 

The most pressing cosmological problem of the SM is perhaps the lack of dark matter 
candidates and another one is a trouble with generating inflation. Both problems may be 
linked to the issue of the instability. Dark matter or an inflaton may come together with 
additional new fields stabilizing the Higgs potential [8-10] and in fact even the Higgs field 
itself may play a nontrivial role in inflationary scenarios [11]. 

The flat spacetime analysis of the stability of the SM is important on its own rights, 
but it may miss new phenomena that arise from the presence of gravity. For example, the 
existence of a non-minimal coupling of scalar fields to gravity which forms the basis of the 
Higgs inflation model [12]. It is worthy to note that such terms are actually needed for 
the renormalization of any scalar held theory in curved spacetime [13, 14]. The problem 
of the influence of gravity on the stability of the Higgs potential was investigated, to some 
extent, using the effective operator approach in [15]. Unfortunately, this approach is based 
on a non-covariant split of the spacetime metric on the Minkowski background and gravi¬ 
ton fluctuations. Its two main problems (apart from the non-covariance) are the limited 
range of energy scales where this split is applicable and the possibility that this method 
underestimates the importance of higher order curvature terms like for example squares of 
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the Ricci scalar, the Ricci tensor and the Riemann tensor. Such terms naturally arise from 
the demand of the renormalisation of quantum field theory in curved spacetime. 

Analyzing Einstein equations with standard assumptions of isotropy and homogeneity 
of spacetime, one can straightforwardly obtain a relation between the second-order curva¬ 
ture scalars (squares of the Riemann and Ricci tensors) and the total energy density. From 
this relation we may see that they become non-negligible at the energy scale of the order of 
10 9 GeV. Therefore, the usual approximation of the Minkowski background metric breaks 
down above such energy. 1 On the other hand, the instability of the SM Higgs effective 
potential appears at the energy scale of the order of 10 10 GeV. This raises a question of 
the possible influence of the classical gravitational field on the Higgs effective potential in 
the instability region. 

Addressing the aforementioned issue is one of the main topics of our paper. To do this 
we calculated the one-loop effective potential for the gauge-less top-Higgs sector of the SM 
on the classical curved spacetime background. We also took into account the presence of an 
additional scalar field that may be considered as a mediator between the SM and the dark 
matter sector. To this end, we used fully covariant methods, namely the background field 
method and the heat kernel approach to calculate the one-loop corrections to the effective 
action. Details of these methods were described in many textbooks, e.g., see [13, 16]. On 
the application side, this approach was used to construct the renormalized stress-energy 
tensor for non-interacting scalar, spinor and vector fields in various black hole spacetimes 
[17-21] and in cosmological one [22]. Recently, it was applied to the investigations of the 
inflaton-curvaton dynamics [23] and the stability of the Higgs potential [24] during the 
inflationary era as well as to the problem of the present-day acceleration of the Universe 
expansion [25, 26]. On the other hand, in the context of our research, it is worthy to point 
out some earlier works concerning the use of the renormalization group equations in the 
construction of the effective action in curved spacetime [27-30]. 

It is important to note that the method we used is based on the local Schwinger-DeWitt 
series representation of the heat kernel (see also [31]), which is valid for large but slowly 
varying fields. In the literature there also exists non-local version of the method engineered 
by Barvinsky, Vilkovsky and Avramidi [32-34] but it is applicable only to small but rapidly 
varying fields. For a more recent development of this branch of the heat kernel method see, 
e.g, [35-37], 

As a final remark we want to point out two other papers that considered the influence of 
gravity on the Higgs effective potential, namely [38] and [39]. In the latter only the tree-level 
potential was considered, while calculations in the former were based on the assumption of 
a flat Minkowski background metric. For this reason, it was impossible there to fully take 
into account the influence of the higher order curvature terms. 

The paper is organized as follows. In section 2 we discuss action functionals for gravity 

1 By definition, for the Minkowski metric we have R = = RaPpvR 1 *^ 111 ' = 0, while for 

the Friedmann-Lemaitre-Robertson-Walker metric R ~ Mp 2 p, R^R^ ~ {Mp 2 p) 2 and R a piivR al3llv ~ 

(. Mp 2 p ) 2 , where p is energy density and Mp ~ 10 18 GeV is the reduced Planck mass. The above 

relations imply that for the energy scale 10 10 GeV we have p ~ (10 10 GeV) 4 and R ~ 10 4 GeV 2 , 

R^R^ ~ Rap^R^ ~ 10 8 GeV 4 . 
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and the matter sector, we also obtain the one-loop effective action in an arbitrary curved 
spacetime. Section 3 is devoted to the problem of the renormalization of our theory, in 
particular we derive the counterterms and beta functions for the matter fields. In section 
4 we ponder the question of the running of the coupling constants and the influence of the 
classical gravitational field on the one-loop effective potential. The last section 5 contains 
the summary of our results. 


2 The model and its one-loop effective action 


As was mentioned in the introduction, in this paper we consider the question of an influence 
of a nontrivial spacetime curvature on the one-loop effective potential in a gauge-less top- 
Higgs sector with an additional scalar field. The general form of the tree-level action for 
this system may be written as 


Sgrav — / y/ i/d X 


1 


(-R — 2A b) + 


167tGb 

+ a± BRiiupaR^ Upa + CH2 bRiwR^ U + OTzbR 2 


( 2 . 1 ) 


S scalar — j[ y/ ffd X (^7 fih B ''j flB ^hB^B^R T ChB^ B hBR ^hB (J^B^'b'J ~f~ 


+ V^XbV^Xb - m 2 XB X 2 B + t XB X 2 B R - \ XB (x) - A hXB X 2 B tf B h B 

*5 fermion / V~gd A x ip B (i^Vg - yBthB^j Vf 


’’B 


, ( 2 . 2 ) 

(2.3) 


where the subscript B indicates bare quantities. 2 

When the scalar interaction term is absent, h represents the radial mode of the SM 
Higgs doublet in the unitary gauge, ip is a top quark and X stands for an additional scalar 
field. The total action is given by 


Stot — Rgrav T Rmat — Sg rav T S sca i ar + Sf erm i on - (2.4) 

To compute the one-loop correction to the effective action we use the heat kernel 
method. Details of the method can be found in [13, 14] (we closely follow the convention 
and notation assumed there). A formal expression for the one-loop correction in the effective 
action is the following: 

f (1) = y lndet (/W 2 D 2 ) , (2.5) 

where \i is an energy scale introduced to make the argument of the logarithm dimensionless. 
In the above relation det means the functional determinant that can be exchanged for the 

2 We used the following sign conventions for the Minkowski metric tensor and the Riemann tensor: 

Vab diag (-\-, , , ), R\rfi d r r T... , R^u R^oa* • 
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functional trace by lndet = Tr In, where Tr stands for the summation over the field indices 
and the integration over spacetime manifold. To find the specific form of the operator D 2 j 
for the one-loop effective action we use the background field method. Our fields have been 
split in the following way: 


X = 


h = 


1 

71 

i 


y/2 ^ 

^ = x + 7 


x + x 

h h 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


where the quantities with a hat are quantum fluctuations and X,h,x are classical back¬ 
ground fields. To find the matrix form of the operator D 2 - we need only the part of the 
tree-level action that is quadratic in quantum fields, namely 


S = 


V^gd 4 x^\X,h,i>}D 2 


X 

h 

$ 


(2.9) 


where we skipped indices of D 2 -. Generally this operator is of the form 


D 2 = □ + U, 


( 2 . 10 ) 


where □ = is the covariant d’Alembert operator and U stands for all non-derivative 

terms. To calculate the one-loop correction we use the following relation [16]: 

T (1) = — lndet(/V 2 T> 2 ) = / V-~gd 4 x tr\^ j — K(x, x, s)|, (2-11) 


where s is a parameter called proper time and K(x,x,s) represents the coincidence limit 
of K(x,x',s). The quantity K(x,x',s ) is the heat kernel of the operator D 2 and obeys 

Q 

i—K(x,x',s) = D 2 K(x,x',s) (2.12) 

os 

with boundary condition lim s _>.o K(x, x', s) = 5(x,x'). In the case at hand, in which fields 
are slowly varying, the heat kernel admits a solution in the form of the Schwinger-DeWitt 
proper time series 


K(x, x', s) = i(47ris) n / 2 exp 


ia(x, x') 
2s 


a vm(x,x') f (x,x',s) 


(2.13) 


where n is the number of spacetime dimensions, cr(x, x') is half of the geodesic distance 
between x and x 1 , A vm is the Van Vleck-Morette determinant 


A vm{x,x') 


-|5(®)l 1/2 |5(^)l 1/2 det 


—d 2 a(x,x') 

dx^dx' v 


(2.14) 
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and F(x, x', s ) = ^2^ = o(is) 3 aj(x, x'), where a,j(x, x') are coefficients given by an appropriate 
set of recurrence relations [16]. Putting all this together the one-loop corrections are given 
by the formula 

r (1) = (2.15) 

where we used the partially summed form of the heat kernel [40, 41] and the trace is 
calculated over the fields (with the correct sign in the case of fermionic fields). The quantity 
fx has the dimension of mass and was introduced to correct the dimension of the action. To 
simplify the notation we introduce M 2 = U + ^R. After integration over s we get 

T (1)=h l ^^ 2(4 l)n/2 ^)a fc M- 2fe }. (2.16) 

Unfortunately, summing the above series is generally impossible. But for our calculation 
we need only its expansion for small s, for two reasons. The first is that beta functions are 
defined by the divergent part of the action which is given by the three lowest coefficients 
(in four dimensions). The second reason is that we are working with the massive slowly 
changing fields for which -C 1. This amounts to discarding terms that are proportional 

to M ~ 2 and higher negative powers of M 2 . Having this in mind we may retain only the 
following terms [40, 41]: 


«0 — lj 

a\ = 0 , 


“ 2 = 1 “ 180 ^^ + Tor\ R i tu P aR ^ Upa + ^ aR \ 1 + * UM2 + 


180 


30 


12 


(2.17) 

(2.18) 

(2.19) 


where W a g = [V Q , Va] (it should be understood as acting on the appropriate component 
of the fluctuation field [X,/i, ?/>]). Using the dimensional regularization we obtain the form 
of the one-loop correction to the effective action 


r (1) ^/« 







+ 2a 2 




( 2 . 20 ) 


where | = 2 — 7 + ln( 47 r), 7 is the Euler constant and n = 4 — e is the number of dimensions. 
Returning to the specific case at hand we find that 


D 2 = 


+ 


-□ 0 0 
0 -□ 

0 0 2 


+ 


~m z x + ixR ~ 3A x X 2 - ^/i 2 
—A hxhX 
0 


—A hxhX 


0 


-m 2 h + ChR ~ 3A h h 2 - ^f-X 2 —jjjVtX 

~75 yth 


' 7 2 VtX 


. ( 2 . 21 ) 
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As one may see, due to the presence of the fermionic field, this operator is not of the form 
□ + U. To remedy this we make the following field redefinitions 3 : 


X —»■ iX, 

h^-ih, ( 2 . 22 ) 

$ n Oreg ] V- 

At this point it is worthy to note that the purpose of the transformation of the fermionic 
variable is to transform the Dirac operator to the second order one. Since the exact form of 
this transformation is arbitrary, it introduces ambiguity in the non-local finite part of the 
effective action as claimed in [42]. This change of variables in the path integral gives the 
Jacobian 


J = sdet 


i 0 0 

0 i 0 

0 0-i (i 7 ^V M 



(2.23) 


where sdet is the Berezinian. For the matrix M that has fermionic (a, j5) and bosonic (a, b) 
entries it is given by 


sdetM = sdet 


a a 
P b 


= det (a — ab P)/det(b) 


(2.24) 


In the case at hand its contribution to the effective action is (omitting irrelevant numerical 
constants) 


J = gl^lndet (*7 M V M -0 reg ) 


(2.25) 


which is proportional to the terms at least quadratic in curvature (R 2 , Ric 2 , Riem 2 ). From 
now on we will work in the limit 6 reg = 0. After the above redefinition of the quantum 
fluctuations the operator D 2 takes the form 



’□ o o’ 


1- 

o 

o 

o 

D 2 = 

0 □ 0 

+ 

0 0 —^ytxn^n 


0 0 □ 


0 0 _ 



+ 


t\R + 3X x X 2 + ^h 2 A hxhX 0 

A hX hX " m 2 h -£ h R + 3\ h h 2 + ^X 2 0 

0 ~72 VtX ~^ R 


(2.26) 


where we used the fact that = □ — ^R. From the relation (2.26) one can see 

that D 2 becomes 


D 2 = U 1 + 2 /i^V m + n, (2.27) 

3 The parameter 8 reg was introduced to ensure invertibility of the considered transformation and it 
should not be identified with the fermionic mass. Moreover, the matter part of the effective action is well 
behaved in the limit of 9 reg —> 0 as can be seen for example in [13]. 
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where 1 is a unit matrix of dimension six and h^ and II are matrices of the same dimension. 
This is not exactly the form of D 2 that we discussed while explaining how to obtain the 
one-loop action via the heat kernel method, nevertheless the formula (2.20) is still valid 
provided we make the following amendments [13]: 


W a p — [V a , V/j] 1 + 2 V[ q /i ( 3] + [h a , hp] 
M 2 = n + ^R1 - - h^. 


(2.28) 

(2.29) 


In the above expression both W and M 2 represent matrices with bosonic and fermionic 
a a 


entries of the form 


P b 


. To take this into account in our expression for the one-loop 


corrections to the effective action we replace tr with str, where str 


a a 

P b 


= tr(a) — tr(b). 


The explicit form of II and hT can be easily read from (2.26), which gives 


M 2 = 


+ 


m 


x 


-{ix-\)R + 2,\ x X 2 + ^h 2 
\ xhX 
0 


\xhX 


0 


m 2 -(Zh-l)R + 3\ h h 2 + ^X 2 0 
0 0 


+ 


0 0 

0 0 

i2y t 




o _ Lr 4- Hl/,2 _ bryV 

V2 x 12^ ' 2 n 2 V 2 » ' . 


(2.30) 


On the other hand, W a p can be computed from the expression (2.28) to be 


W a p 


0 0 0 

0 0 ~2V2 - V/3X7a) - yhx (7a7/j “ 7/37a) 

0 0 \R a pyuYY + ^ - Vphxa) - \h 2 (7 q7/3 - 7s7a) 


(2.31) 


To summarize the calculations, we present below the full form of the renormalized one-loop 
effective action (T) for the matter fields propagating on the background of the classical 
curved spacetime. The details of the renormalization procedure will be given in the next 
section. In agreement with our approximation, we keep only the terms proportional to the 
Ricci scalar, its logarithms, the Kretschmann scalar {R a pii^R al3>lu ) an d the square of the 
Ricci tensor. Moreover, we discard terms proportional to the inverse powers of the mass 
matrix and renormalize the constants in front of higher order terms in the gravity sector 
to be equal to zero (07 = 07 = 03 = 0) at the energy scale equal to the top quark mass. 
Additionally, we disregard their running since it is unimportant from the perspective of the 
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effective action of the matter fields. The final result is 


r = — 


1 


I67 tG 


j y/^gd*x(R + 2A) + J yf^gd A x jx 


- -^yh 


X+ 


+ £ h R) h 2 - ^ ^ h fh 2 X 2 + 

+ l -V^X - i (nix ~ txR) X 2 - ^X 4 + 


+ 


647T 2 


VtX (i'fVp + 2-^=y t hj x ~ ^y / 2 V /J ltV , '7i - 2 \y 2 t In 


X p hX 7 v li+ 


l l 


+ 


+ 


- |fr(o>to (4j)) + fain (4,) - tr (a 2 In (4;)) + f *™ 2 + 8i> 2 l n ( 

+ ln {~^) R ~ y t h4 ln 

— 180 (~ Ra P Ral3 R apiiv 

-\R a ^R a ^ ln(^)]} 

where a and b are given by 

b = ^< 2 '* 2 - h R - 




- 12 b 2 + 


ln, ^) +ln (^)- 21n (^ 


a = 


m% 


-{ix-\)R + 2,\ x X 2 + ^h 2 


\xhX 

The eigenvalues of the matrix a are 


m, 


^hxhX 

-(Z h -l)R + 3\ h h 2 + ^X 2 


1 

fl± = 2 


m x + m h ~~ ( €x + £h ~ g j R + f 3A^ + — Xhx J h 2 + ^3 Ax + -A hx ) X 2 


(2.32) 

(2.33) 

(2.34) 

+ 


m 2 x — mi—[t,x — £h)R+[ x^hx ~ 3A h )h 2 + [ 3Xx — -A hx X 2 


1 


+ 4 (a hX hx' 


(2.35) 


3 Divergent parts of the one-loop effective action and beta functions 
3.1 Divergences in the one-loop effective action 

Divergent parts of the one-loop effective action of our theory can be straightforwardly read 
from the expression (2.20) and they are given by the sum of terms proportional to -. In 


our case 


di) _ 

div 


= / y/^dTx- 


4 i n 2 


£ (47r) 2 4 

+^OM 2 + ^W af} W ay | 


- < strM A + 2 sir 


- RuvR* aU 4 - Riivorj 

180 p 180 p p 


Riivp ° + ^o nR * 1+ 


(3.1) 


The terms proportional to □ are full four divergences and can be discarded due to the 
boundary conditions. Moreover, we also neglect the terms that are of the second and 



















higher orders in the curvature, since they contribute only to the renormalization of the 
gravity sector. A precise form of this contribution is well known and can be found for 
example in [14]. Having this in mind the only relevant terms are strM 4 and strW 2 . We 
may write the matrix M 2 in the following form: 


M 2 


a 



+ l 4 h x 


(3.2) 


where a and b were defined in (2.33) and (2.34). Having this in mind, the only relevant 
entries of M 4 are the diagonal ones (off-diagonal entries do not contribute to sir) 


M 4 




a 

v hY - i^=y t bV' 


(3.3) 


From this we obtain 


strM 4 = tr(a 2 ) — y 2 \ 


- 2 ~^yth 


X~8b 2 + yiV.hVh. 


(3.4) 


In the above expression we have doubled fermionic contributions to restore proper numerical 
factors changed due to nonstandard form of the fermionic gaussian integral used by us. The 
second term that contributes to the divergent part of the one-loop effective action comes 
from 


2 str 2-W a yW a ^ = pstr( 
12 D 


0 a 

0 c 2 


(3.5) 

(3.6) 


with 

c 2 = ^Ra^R aP pa'f'f'f'f - ^y 2 t h 2 ( 7a7/3 - 7/37q ) R ap 

- ly 2 (V«/l 7j9 - Vfih'Ta) (v“^ - V^ 7 “) - 2-y 2 h 2 R Q ^Y ( 7 V - 7 V) + 

+ ^yt h * ( 7 « 7 /3 - 7 /3 7 a) ( 7 “ 7 /3 - W) , (3-7) 

where we omitted the terms proportional to the odd number of gamma matrices. Combining 
above expressions with the one for the M 4 divergent part of the one-loop effective action 
gives (after discarding purely gravitational terms of the order 0(1Z 2 ) , where 7 Z 2 represents 
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terms quadratic in Ricci scalar, Ricci and Riemann tensors) 


divp r (1) = 
divp r (1) = 


2 h 


e 647r 2 ( 

2 h 


+ 


e 647r 2 
2 


str(M ) - -tr(c) j , 

m\ - (£ x - + 2>\xX 2 + 

6 2 


(3.8) 


+ 


1 , 


Aft,X 


m h — (£/i — q)R + 3Aft/l + ^ 


X 2 


+ 2 rfixh 2 X 2 - y 2 t x -i^X v - 2^=y t h 


1 

'3* 


+ 

X -2 


y 2 h 2 — -R 

yt 6 


+ Vth 4 + 


1 


+ 2yfX l _ l hX IJ h-Xy 2 h 2 R}. 


(3.9) 


3.2 Counterterms and beta functions 

After finding the divergent part of the one-loop effective action we shall discuss the renor¬ 
malization procedure in detail. The matter part of the tree-level Lagrangian in the terms 
of bare fields and couplings can be written as 


Cb matt = \v a h B V a h B - i [m 2 hB - £ hB R] 


h%- 


A hXB 


h 2 B X 2 B + 


+ l -X a X B V*X B - i \m 2 XB - £ XB R\ X 2 b - ^fx%+ 


+ XB 


h B 


XB- 


(3.10) 


The same Lagrangian can be rewritten in terms of renormalized fields and coupling con¬ 
stants. Appropriate relations between bare and renormalized quantities are 

h B = zl /2 hfT*, X B = Z]l 2 Xpr §, xb = z][\p-\, 

A hB = Z h 2 Z\ h p £ Xh, A xb = Z x Z\ x p £ X\, XhXB = Z h 1 Z x Z\ hx p £ Xhx, 

VtB = Z~ 1 Z~ l/2 Z yi Z2 £ y, 

m hB ~ Z h Z mh m h , m XB — Z x Z mx rn x , 

ihB = Z^Z^ h , £xb = Zx'Zfrtx, ( 3 - 11 ) 

where we introduced mass scale p to keep quartic and Yukawa constants dimensionless. 
Using the above formulae and splitting the scaling factors as 2 a = 1 + 5 a we may ab¬ 
sorb divergent parts of one-loop corrections to the effective action divp T^. One-loop 
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counterterms are 


1 h 


SZ h = - 

i h 


£ (4ir) 2 


= ~ 

SZx x = 
bm\ = - 

^h = 

<5£x = 


£ (47t) 2 
l h 

£ (47r) 2 
l h 

£ (47r) 2 
l h l 

£ (4vr) 2 £ h 

l h l 


2 2/t, 5Z X = 0, 

1 \ 2 'll 4 

l 8Ah + i^-2^ 


18A X + 


2 A h A h 

1A? 


XV _ 1 ^ 4 2 

5 X ~ £ (47t) 2 2^ t ' 

sz„ = - h 


'hX 


2 Ax 


^ = 


V £ (47t) 

1 h 


2 yl 


£ (4ir) 2 


[6Aft + 6Aj; + 4 \ hX ], 


6A h + A hx —f 


m 


h - 


5m 2 x = - h 


£ (47t) 2 

1' 


6 Ax + A/,x- 


mi 


m 




£ (4vr) 2 £x 


6Afe ( & — g J + A/jx ( £x - 6 

6Ax Ux - + Aftx Uh - ^ 


;Vt 


(3.12) 


Using the above form of counterterms we may compute beta functions for the quartic and 
Yukawa couplings 


o ^ 5 o 

Pyt / i_\9 nVt • 


yt (4vr) 2 2 

h 


h h = 

P\x = 


h (4vr) 2 

h 


P\ hx = 


(47t) 2 

h 


18A| - 2 yf + 4 y 2 X h + -A| x 


18Ax + ~2^hx 


[4A^x + 6A/ l x(A/ l + Ax) + 2A hxy 2 


AhX (4vr) 2 

Analogous calculations give us beta functions for masses and non-minimal couplings 


f> 2 =^- 

Pm l (4vr) 2 

h 


(4ir) 2 

h 

(47r) 2 

h 

(47r) 2 


[6A h m 2 h + 2 y 2 m 2 h + A ft xm|] , 


P m \ = 7 X ^2 [ 6X ^ m x + A hxm\ 

hn = 

hx = 


6A/i(6i — -) + X hX (£ x — -) + 2 y 2 {ih ~ x) 
bob 


6Ax(£x - x) + Ahx(?/i - x) 
b b 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


For completeness, we also give the anomalous dimensions for the fields (computed according 
to the formula 70 = ^^f) 


h 


ih = 


(471< 

7x = 0, 


2 yt> 


7x = 


h l 


(47t) 2 4 


Vt- 


(3.21) 

(3.22) 

(3.23) 
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At this point we can compare our results for the beta functions of the nonminimal couplings 
of the scalars to gravity with those obtained for the pure Standard Model case [24]. If 
we disregard the modification of & stemming from the presence of the second scalar, 
namely the A hx component, we are in agreement (modulo numerical factor due to different 
normalizations of the fields and the absence of vector bosons in our case) with the results 
from the cited paper. 


4 Running of couplings and stability of the effective scalar potential 
4.1 Tree-level potential and the running of the couplings 


Our theory consists of two real scalar fields (corresponding to the radial mode of the Higgs 
scalar in the unitary gauge and an additional scalar singlet) and one Dirac type fermionic 
field that represents the top quark. From now on, we will call the second scalar the (heavy) 
mediator. To solve the RGE equations for our theory we need boundary conditions. A 
scalar extension of the Standard Model was extensively analyzed in the context of recent 
LHC data (for up to date review see [10]). We use this paper to obtain initial conditions 
for RGEs of the scalar sector of our theory. An energy scale at which these conditions were 
applied has been set to /it = 173 GeV. 

Form of the tree-level potential 


V T ree(h, X) = i (m 2 x - £ X R) X 2 + ^A 4 -I- ^h 2 X 2 + 
+ 2 ( m h - £hR) h 2 -I- -j-h 4 


we may find the tree-level mass matrix 


M 2 = 


m 


x 


- ixR + 3A X A 2 + ^h 2 

\xhX 


\xhX 
m 2 -C h R + 3X h h 2 + ^X 2 


(4.1) 


(4.2) 


At the reference energy scale /it the Vr r ee{h, X) has one local maximum h = 0, X =0, two 
saddle points h = 0, X ^ 0 and h ^ 0, X = 0 and one local minimum 

Xh.x 


m h — €hR + H X~ — 0, 


m x — £,xR- + XxX + 


2 , ^hX u 2 


/r = 0. 


(4.3) 


We identify this minimum with the electroweak minimum (electroweak vacuum) where the 
mass matrix (4.2) takes the form 


M 2 


2X x X 2 X hX hX 

XhxhX 2A hh 2 


(4.4) 


Replacing the fields by their physical expectation values h = Vh,X = v x we may define 
physical masses as the eigenvalues of the above matrix 

m 2 H_ = ^ X x v x + X h v 2 h - \J {X x v 2 x - X h v 2 h ) 2 + (A^y^.y^) 2 ^ , (4.5) 

m 2 H+ = ^X x v x + X h v\ + \J ( X x v 2 x - X h v 2 h ) 2 + (Xhxvxvh) 2 ^ . (4.6) 
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For concreteness, we set these masses to mu_ = 125.5 GeV and mn + = 625 GeV. This 
choice amounts to identifying the lighter of the mass eigenstates with the physical Higgs 
and the heavier one with the scalar mediator outside the experimentally forbidden window. 
Moreover, we take vev of the Higgs to be Vh = 246.2 GeV. The expectation value of the 
second field can be expressed by the parameter tan(/3) = This parameter is constrained 
by the LHC data to tan(/3) < 0.33 for the m# + < 700 GeV [10], we fix it to the value 
tan(/3) = 0.33. From the Lagrangian of the scalar sector of the theory at hand one can 
see that it is described by five parameters, namely two masses and three quartic couplings. 
So far, we specified four parameters: two masses (mass eigenstates) and two vevs of the 
scalars so we have one more free parameter. This parameter is the mixing angle between 
mass eigenstates i7_ and H + and the gauge eigenstates X and h 



cos a — sin a 
sin a cos a 



(4.7) 


We fix it as sin(a) = 0.15. Remembering that the above rotation matrix diagonalizes the 
mass matrix (4.4), we find an explicit expression for the mixing angle a (— ^ < a < -|) 


sin 2a = 


cos 2a = 


_A hxVhVx _ 

\J (A hvl - Xxvx) 2 + {Xhxvxvh) 2 
_ X xv 2 x - A hvj _ 

yj (A hv\ - \xv\ ) 2 + (A hxvxvh) 2 


(4.8) 

(4.9) 


Using the formula for the sin(2a) and relations (4.5) we express the quartic couplings in 
terms of physical masses, vevs and a 


A h = 
Ax = 
A hx = 


m 


H- 


m 


H 


+ 


H a 


— m 


H- 


H 


sin 2 (a) 


nrii 


m. 


m H _ m H + - m H _ 2 / \ 

H-o-cos (a) 


2v x 


2v x 


H 


m f 

2v x 


■ cos 2 (a) H-j- sin 2 (a) 


H 


sin 2 (a) + 


m. 


2v x 


cos 2 (a) 


m H+ - m H _ 

2v h v x 


sin(2a) 


(4.10) 

(4.11) 

(4.12) 


In figure 1 we present the dependence of the quartic coupling on the mixing angle for 
tan(/?) = 0.33 and physical masses and vevs chosen as stated above. From this plot we may 
infer that as we increase the mixing angle the Higgs quartic coupling (A h) and the interaction 
quartic coupling (A hx) increase while the heavy mediator quartic coupling (Ax) decreases. 
Moreover, there are two non-interacting regimes. The a = 0 case corresponds to the 
two non-interactive scalars with quartic selfinteraction. Additionally, from the form of the 
beta function for A/jx we may infer that the interaction between these two fields will not be 
generated by the quantum corrections at the one-loop level. The second regime corresponds 
to a = |, but since for this value of the mixing angle also Ax is zero the additional scalar 
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is tachyonic. As the final step we express mass parameters of the Lagrangian in terms of 
our physical parameters. To this end, we use equations (4.3) and obtain 


~ m x ~ ^xv'x 4 2 ~ v h ~ (4-13) 

~ rn h = ^hv\ H —v\ ~ ChR- (4-14) 


The remaining parameters of the scalar sector are the values of the non-minimal coupling 
to gravity ^ and £x and the field strength renormalization factor for the h field. We have 
chosen Zh = 1 at the reference energy scale yt and for a nonminimal coupling we considered 
two different cases. The first one was Ch = Cx = 0 which results in and £x becoming 
negative at high energy. The second one was £/i = £ x = | for which and £x stay positive 
at high energy. The choice of the initial conditions was arbitrary but allowed us to present 
two types of the behavior of the running of the nonminimal couplings, that will be discussed 
shortly. 

Having fully specified the scalar sector, we now turn to the fermionic one. It possesses 
two parameters, namely the field strength renormalization factor and the Yukawa coupling 
constant. The first one is naturally set to unity at yt and we set the top Yukawa coupling 
as yt = 0.9359, where the physical top mass was chosen as mt = 173 GeV. In figure 2 
we present the running of the Yukawa and quartic couplings. For the described choice 
of the parameters point, at which the Higgs quartic coupling becomes negative, is given 
by t « 52.2, which corresponds to the energy scale y ~ 10 10 GeV. We also observe that 
the most singular evolution will be that of the heavy mediator quartic coupling Ax and, 
indeed, this coupling hits its Landau pole around t ~ 59. In figure 3 we depicted the 
running of the mass parameters of the scalar fields. This running is quite big and amounts 
to an increase of more than 50% around t ~ 50. Figure 4 presents the running of the 
non-minimal couplings and field strength renormalization factors. In figure 4b we may see 
that the non-minimal couplings run very mildly and they are positive for the low-energy 
region and become negative for high energy regions. Just to remind, our initial condition 
for them was £/, = £x = 0 at the reference point y = mt = 173 GeV. On the other hand, if 
we choose initial values for ^ and £x to lay above the so-called conformal point £ = g, they 
will stay positive in the whole energy region. This type of behavior is visible in figure 4c. 


4.2 One-loop effective potential for scalars in curved background 

In this subsection we present the form of the one-loop effective potential for the Higgs- 
top-heavy mediator system propagating on curved spacetime. In the framework of the 
R-summed form of the series representation of the heat kernel (the subset of the terms 
proportional to the Ricci scalar is summed up exactly) and on the level of the approximation 
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Figure 1: The dependence of the initial value of the quartic couplings on the mixing angle 
a. 
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Figure 2: The evolution of the Yukawa and quartic couplings of the scalar fields. The 
running scale range is from //q — 2.7 K to n max = 10 11 GeV. 
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Figure 3: The running of mass parameters for the scalar fields. The energy range is from 
Mo = 2.7K to n m a X = 10 11 GeV. 



- m h 2 lm h 2 (ij t ) 

. m x 2 /m x 2 {iJt) 


discussed earlier we may write it as 

Vone-ioop = - j - \ {ml - ChR) h 2 - ^h 4 - h 2 X 2 - 1 (m 2 x - txR) X 2 - ^x 4 + 



where b is given by (2.33) and a± is defined by (2.35). Let us recall that in our approximation 
we are discarding terms of the order 0 (^-), where a.j = {a+,a_, 6 } and 1Z 3 stands for all 
possible terms that are of a third order in curvature. Since we specialize our considerations 
to the cosmological case we take the background metric to be of Friedmann-Lemaitre- 
Robertson-Walker type, for which 


R = —6 



(4.16) 


where A is the scale factor, namely the metric is ds 2 = dt 2 — A 2 (dx 2 +dy‘ 2 +dz 2 ). Meanwhile, 
the Einstein equations reduce to the so-called Friedman equations 


2) 

/ 

(4.17) 

2 j + H 2 = -Mp~ 2 p , 

(4.18) 
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(a) The running of the field renormalization factors for h and the 
fermion field X- The energy range is from fj-o = 2.7 K to fj-max = 
10 11 GeV. 



(b) The running of the non-minimal couplings to the gravity for the 
scalar fields, the initial conditions were £/, = £x = 0 at the /r = mt- 
The energy range is from fi o = 2.7 K to Umax = 10 11 GeV. 



(c) The running of the non-minimal couplings to the gravity for the 
scalar fields, the initial conditions were £h = £x = | at the [i = mt- 
The energy range is from fj-o = 2.7 K to [imax = 10 11 GeV. 

Figure 4 
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— _2 

where Mp = 8nG is the reduced Planck mass, p is energy density and p is pressure. 
Using the above equations we may tie the Ricci scalar to the energy density and pressure 


R = -3 M P 2 


-p+\p 


(4.19) 


The other useful scalars are 


-R a pR aP + R a ^R a ^ u = ~ 12h2 J = 2 m p A P Qp + p) 


= 3 \Mp 2 P ) , 


R a p^R a ^ v = 12 


" 4 +'iy 


= 12 M 


-4 


1 o 1/1 


9 P + 4[3 P + P 


= 3(^P 2 P) , 


(4.20) 

2 


(4.21) 


where the last equalities are valid in the radiation dominated era, where p = \p- Now let 
us discuss what the above statements mean in the context of our approximation. From 
relations (4.20) we may infer that 7£ 3 ~ (Mp 2 p) 3 . On the other hand, our expression 
for the one-loop effective action is valid when terms that are of higher order in curvature 
are suppressed by the held dependent masses a+, a_ ,b. This means that in order to 
investigate the influence of a strong gravitational held on the electroweak minimum (small 
helds region) we must have ^ <C 5,2, where m ~ 10 2 GeV is the mass scale. This leads 
us to the relation ~ 10 2 4- 10 3 GeV 2 . To connect the energy density to the energy 
scale we use the formula p = ow 4 + /i 4 , where p is the running energy scale (as introduced 
in RGE) and a is a numerical constant. For our particular choice of p = -^=/i we have 
p = ctp 4 + (^^) 4 - We choose v = 10 9 GeV and a in such a way that our approximation is 
valid at the electroweak minimum. 

At this point an additional comment concerning the running energy scale is in order. 
In the previous section we presented results concerning the running of various couplings 
in the model. To this end we considered the energy scale present in RGEs as an external 
parameter, but this is sometimes inconvenient for the purpose of presenting the effective 
potential. For this reason, in this section we adopt the standard convention (in the context 
of studies of the stability of the Standard Model vacuum) of connecting the energy scale 
with a held dependent mass. In the theory at hand this leads to the problem of a non¬ 
uniqueness of such a choice since we have three different mass scales (a+, a_, b). To make 
our choice less arbitrary we follow some physical guiding principles. First of all, the energy 
scale should be always positive. Secondly, the relation between helds and the energy scale 
should be a monotonically increasing function. This condition ensures that an increase of 
the value of helds leads to an increase of the energy scale. Moreover, we expect that for 
a single given helds configuration we get a single value of the energy scale. Having this in 
mind we discard a+ and a_, because they are not monotonic functions of the helds. This 
leads us to the choice p = b = -^h, where we discard the gravity dependent term R since 
it is zero at the radiation dominated era. 

After explaining the choice of the running energy scale in more detail, we want to elab¬ 
orate on the physical meaning of the connection between the total energy density and the 
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running energy scale. At the electroweak minimum we still may observe large gravitational 
terms due to the fact that most of the energy is stored in a degree of freedom other than 
the Higgs field, this is represented by the constant (field independent) term v. On the other 
hand, we expect that in the large field region h> u a significant portion of the total energy 
density will be stored in the Higgs field itself (in the scalar field sector in general). The 
amount of this portion is controlled by the parameter a. 

Since the reduced Planck mass is of the order of 10 18 GeV, this leads us to the conclu¬ 
sion that the maximum energy scale at which our approximation to the one-loop effective 
potential around the electroweak minimum is valid is of the order u 10 9 GeV. Above this 
energy scale terms that are of higher order in curvatures become large and we need another 
resummation scheme for the heat kernel representation of the one-loop effective action. It 
is worthy to stress that despite the fact that the finite part of the one-loop action becomes 
inaccurate above the aforementioned energy scale, the running of the coupling constants is 
still described by the calculated beta functions. This is due to the fact that the UV diver¬ 
gent parts get contributions only from the lowest order terms in the series representation 
of the heat kernel. 

We plot the one-loop effective potential in the Friedmann-Lemaitre-Robertson-Walker 
background spacetime in figure 5. The total energy density that defines curvature terms 
was set as p = aix 4 + /i 4 , where a = 50, v = 10 9 GeV and [X = Figure 5a represents 

the small field region (the region around the electroweak minimum). For given parameter 
choices the expectation values of the field are Vh = 246.2 GeV and vx ~ 746 GeV. The 
black line in this figure represents the set of points in the ( X , h ) plane for which our one- 
loop approximation breaks down. For these points one of the eigenvalues of the matrix (3.2) 
becomes null and terms (discarded in our approximation as subleading ones) proportional 
to the inverse powers of this matrix become singular. Moreover, to the left of this line the 
one-loop potential develops an imaginary part due to the presence of the logarithmic terms. 
In figure 6 we present the influence of the gravity induced terms on the effective potential 
in the radiation dominated era. To make the aforementioned influence of the gravitational 
terms clearly visible, we choose a single point in the held space. Namely, we choose the 
electroweak minimum for which h = Vh and X = vx ■ We may see that for large total energy 
density this minimum becomes shallower. In figures 7 and 8 we also plot this influence for 
other cosmological eras, namely matter dominated and the de Sitter ones. From figure 8 we 
may infer that for the de Sitter era and positive £ the minimum becomes even more shallow 
and this effect is orders of magnitude bigger than for the radiation dominated era. This is 
mainly due to the fact that terms which contribute most are from the tree-level part of the 
effective potential. On the other hand, from figure 7 we infer that for the de Sitter era and 
£ = 0 the one-loop gravitational terms (tree-level ones are zero due to £ = 0 ) lead to the 
deepening of the electroweak minimum. The magnitude of this effect depends on the total 
energy density. 

Another interesting question is how big should the gravity induced parts be to qualita¬ 
tively change the shape of the effective potential. To get the order of magnitude estimate 
we consider only the Higgs part of the effective potential. For now, we specify the back¬ 
ground to be that of the radiation dominated epoch, for which we have p = \p and R = 0. 
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(b) 

Figure 5: The one-loop potential for the scalar fields in (a) small and (b) large field 
regimes. The running energy scale was chosen as /i = cr = 50 and u = 10 9 GeV. The 
thick black line in (a) represents a set of points for which a_ = 0. The dashed line in (b) 
represents the line along which V 1 = 0. 
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Figure 6: The influence of the gravity induced terms on the one-loop potential for fixed 
values of fields (h = Vh,X = vx)- The running energy scale was set to be equal to the top 
mass p = mtop■ Total energy density is given by p = <rpf otal = av A + mf op . 
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Figure 7: The influence of the large curvature on the electroweak minimum for various 
equations of state: rad - radiation dominance (p = ^p), dS - de Sitter like (p = —p), matt 
- matter dominance (p = 0). The energy density was given by p = ow 4 + yU 4 , where a = 1 
and p = The non-minimal couplings were £/i = £x = 0 at p = mt- The insert shows 
a close up of the behavior of the gravitational corrections for the radiation dominated era. 
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Figure 8: The influence of the large curvature on the electroweak minimum for various 
equations of state: rad - radiation dominance (p = ^p ), dS - de Sitter like ( p = —p), matt 
- matter dominance (p = 0). The energy density was given by p = a(v 4 + p 4 ), where a = 1 
and p = The non-minimal couplings were = E,x = i at p = rrq. The insert shows 
a close up of the behavior of the gravitational corrections for the radiation dominated era. 



rad 

dS 

matt 


In the small field region the most important fact defining the shape of the potential is the 
negativity of the mass square term m 2 < 0. Meanwhile, gravity contributes to the following 
terms: 


I/W = 1 

9rav 64vr 2 


4 

180 


-R af) R aP + Ra8auR a ^‘ 


'OlPlM/-* 


i „(^)+to( 


p 2 


- 2 In 




+1^^^ (A 


+ 

(4.22) 


With our convention for the running energy scale we see that the fermionic logarithm 
In is equal to zero and the remaining two logarithmic terms are positive and of the 

order of unity, see figure 9a. Let us call their total contribution b. Having this in mind, we 
may write 


V{h 2 ) = 


\m 2 h h 2 + 1 


; m h + 


64tt 2 180 
1 4 

64vr 2 180 


-RcfiK* + Ro:0^R a>6lIU ) b = 
(~R a pR aP + R a ^R a ^) ■ 


h 2 = 


—m h + 
2 


1 4 4 


64vr 2 180 3 l Mp2p j h 2 


, 2 2 7,2 

h = m e ffh , 


(4.23) 




where we put h = 1. Since we are interested in the influence of the gravity on the electroweak 
minimum, we make the following replacement in the bracket: h = Vh, also for the chosen 
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physical Higgs mass, vev and mixing angle we have = —6.1 • 10 4 GeV 2 . Now our goal is 
to determine the energy density at which m 2 e ff >0. It is given by 


135 

p = 4xv h \m h \\j Mp 


(4.24) 


and corresponds roughly to the energy scale v ~ 10 10 4-10 11 GeV (under the assumption of 
p = z/ 4 ). This value is slightly above the energy scale for which our approximation is valid 
(in the case of the small field region), nevertheless it is reasonably below the Planck scale. 
For the de Sitter case, on the other hand, the dominant contribution comes from the tree- 
level term representing the non-minimal coupling of the scalar to gravity. Straightforward 
calculation gives 


P = ^ M p\m 2 h l (4.25) 

which leads to the energy scale of the order of v 10 10 GeV. 

Before we proceed to the large field region we want to consider the temperature de¬ 
pendent correction to the effective potential. Specifically, we will focus on the influence 
of the curvature induced term on the critical temperature for the Higgs sector of our the¬ 
ory. The leading order temperature dependent terms in the potential will contribute as 
htemp ~ $h 2 T 2 , where (3 is a constant that depends on the matter content of the theory. 

First, let us focus on the beginning of the de Sitter era when most of the energy is 
still stored in the fields excitations. 4 In this case we may assume T = v, p = <rzH and the 
dominant curvature contribution comes from the tree-level non-minimal coupling term 


V(h 2 ) = 


-mi+(3T 2 --£, h R 


h 2 = 


+2 £ h M p au 


(4.26) 


where we used the Einstein equation to express the Ricci scalar by the energy density 
and assumed that ^ > 0. From the above relation we may find the critical temperature 
(critical energy scale i/ c ), for which the origin becomes stable in the direction of h. After 
some algebraic manipulation we obtain 


= 5 Mp 


P 


+ 


P 


2 2£/i<t 


'1 + 


4 £fetT|m,fcl 

P 2 Mp 


(4.27) 


Expanding the square root in the last equation in its Taylor series and relabeling v c = T c , 
we get the following formula for the critical temperature: 



mi\ 2 £ h (T 


(4.28) 


P 2/3 3 Mp 

where we keep only the first two terms in to the Taylor series. Comparing it with the flat 
spacetime result T£ lat = y we can see that the gravity contribution is suppressed 

4 We did not consider the preinflationary era but the short de Sitter period in the middle of the radiation 

dominated era that sometimes is introduced to dilute the relic density of the dark matter. 
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by Mp 2 . Actually, this result also applies to the matter dominated era (up to a numerical 
factor that stems from the modification of the relation between R and p in matter dominated 
era). 

On the other hand, deep in the de Sitter era the energy stored in matter fields is diluted 
by the expansion and the only relevant source of temperature is the de Sitter space itself 5 
(this amounts to setting /3 = 3A/, in first part of (4.26)). This temperature is given by 
T dS = which can be expressed through the Einstein equations by the energy density as 

T dS = Using the last relation to express the energy density by temperature and 

2v37t 

plugging the result into (4.26) one obtains 


ndS 



1 

\X h + 87t 2 £/i| 


(4.29) 


This expression gives the critical temperature above which the electroweak minimum be¬ 
comes unstable. It is interesting to note that, contrary to the previous case, the gravity 
contribution is multiplicative and inversely proportional to the non-minimal coupling con¬ 
stant £/j. This implies that if ^ is big, like for example in the case of the Higgs inflation 
where it is of the order of 10 4 , the critical temperature may be an order of magnitude smaller 
in comparison to the one calculated with the assumption of flat background spacetime. 

As the next case we consider the radiation dominated era. To find the critical temper¬ 
ature we need to solve the equation 

+ 6i^I5o {- R ^ Ral> + ^^ + 4r 2 = o, (4.30) 


where b is defined as in (4.22). Using Einstein equations to eliminate the squares of the 
Riemann and Ricci tensors, assuming T = v, introducing a new variable x = u 2 and defining 


a small coefficient ckq = jfo we may rewrite the above equation as 


r— 4 


a 0 x 4 + fix - - \m 2 h \ = 0. 


(4.31) 


The formulae for the general roots of the fourth order polynomial are quite unwieldy and 
can be found for example in [43]. Using Mathematica computer algebra system we found 
that this equation possesses only one real positive solution, with a series representation (the 
Maclaurin series in qq) given by 


mi 


x 


2/3 16/3 5 


\m 2 h \ A , rM 5/34 

—4- a 0 + O(« 0 ) 


(4.32) 


From the above relation we find the critical temperature for the radiation dominated era 


T — 

± C — 


mi 


2/3 


\m 


214 


1 16& „■> 
M 


-4 


16/3 5 647T 2 640u^ 


(4.33) 


5 The Hawking temperature T dS enters through de Sitter fluctuations of the scalar field substituted into 
the quartic term in the Higgs effective potential. 
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The first observation is that the gravitational terms induce only an additive correction 
to the critical temperature. The second one is that this correction is suppressed by the 
factor Mp 4 so its influence on the aforementioned temperature is very small. This is in 
contrast with the de Sitter case where the gravitational correction may, in principle, change 
the temperature even by an order of magnitude due to the multiplicative nature of these 
corrections. 

Now we turn our attention to the large held region. The most important term of the 
potential is : ^p-/i 4 , where \ e ff contains factors coming from the running of the Higgs quartic 
coupling and the usual held dependent parts coming from the one-loop correction (in the 
absence of gravity). Taking the gravity into account, the relevant part of the potential is 


V{ tf) = hiM h 4 + hW„. (4.34) 

From figure 9b we may see that in the large held region h ~ 3 ~ 4 • 10 10 GeV all logarithms 
are of the order of unity. Although all logarithms are roughly of the same order, the 
leading contribution comes from the fermionic one. This is due to the fact that in Vg}\ v 
the contributions dependent on a± are multiplied by the prefactor that is ten times smaller 
than the term ^R a /3^ u R a ^^ 1 ' hi Now we may write the relevant part of the effective 

potential 


V(^) = + (£) = 


4 4; 


^eff(h) + 64^2 3 3 \M p 2 p 


c 

h 4 


- \h=ho 


h 4 = \x eff {h)h\ 


(4.35) 


where c = In is a number of the order of unity. In the large held region ho ~ 3T0 10 GeV 

we expect that A e //(ho) = d < 0. Now we want to address the issue of how big should the 
energy density be in order to make X e ff(ho ) positive. The straightforward calculation gives 


P = 


AirhyMp 



(4.36) 


For d = \X e ff\ r^j 0.02 we obtain the energy scale v 10 14 GeV. This is again slightly 
above the region of validity of our approximation (which is v rs j 10 10 GeV for the large 
held regime), but still much below the Planck scale. Turning again to the de Sitter era we 
find that the dominant contribution comes from the non-minimal coupling of the scalar to 
gravity. Writing the relevant piece of the potential as 


V(h‘) = T 


x eff(h) - 


h\ 


\h=ho 


we may deduce that the critical energy density is given by 

P = ^Ihld. 


(4.37) 


(4.38) 
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Figure 10: The effective quartic Higgs coupling, as defined by the relation A h e ff(h) = 
; for various equations of state: flat - flat spacetime result, rad - radiation 
dominance ( p = -|p), dS - de Sitter like ( p = — p). The energy density was given by 
P = Phc + (^7§) 4 > where Phc was specified by the relation (4.36) and equal to ph c = 
(2.04 • 10 14 GeF) 4 . The X field was constant and set as equal to X = vx- The non- 
minimal couplings were ^ = 0 at the p = rnt . The insert shows a close up of the 

difference between the flat spacetime and the radiation dominated era. 


In the above formula ho and d are defined in the same manner as for the radiation dominated 
era. For the same value of ho and d like in the previous case we obtained the following 
energy scale at which the discussed effects are important: u 7 ■ 10 13 GeV. Obviously, if 
£h becomes negative, for example due to the running (figure 4b), we always get worsening 
of the stability, \ e ff becomes negative for the lower energy scale than in the flat spacetime 
case. The discussed effects are illustrated in figures 10 (for negative £) and 11 (for positive 
£). Although the obtained energy scales seem to be high (for both radiation dominated 
and the de Sitter eras), the associated energy density is of the order p ~ ICG 21 4- 10 ~ 20 pp, 
where pp = Mp is the Planck energy density. 

Figure 5b presents the large field region of the effective potential. The thick dashed 
line represents a set of points for which V = 0. Below and to the right of this line the 
effective potential becomes negative which indicates the region of instability in the field 
space. This region starts around the point (X = 0G eV,h ~ 4 • 10 10 GeV) and expands 
towards the larger values of h and X fields. In figure 12 we depicted the effective potential 
one-dimensional trajectory in the field space starting at the electroweak minimum and 
ending in an instability region. For this purpose we fixed values of X field by the following 
conditions: X = vx or X = Ah + B. In the latter case the coefficients A and B were chosen 
in such a way that the straight line connects points (vh,vx) and (h m , 0), where h m lies in 
the instability region. From the discussed figure we may infer that the actual trajectory 
connecting the electroweak minimum and the instability region is not very important. The 
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Figure 11: The effective quartic Higgs coupling, as defined by the relation A h e ff(h) = 
, for various equations of state: flat - flat spacetime result, rad - radiation 
dominance (p = -|p), dS - de Sitter like ( p = — p). The energy density was given by 
P = Phc + (^7§) 4 > where Phc was specified by the relation (4.36) and equal to ph c = 
(2.04 • 10 14 GeF) 4 . The X field was constant and set as equal to X = vx- The non- 
minimal couplings were £/j = fx = | at the p = mt . The insert shows a close up of the 
difference between the flat spacetime and the radiation dominated era. 
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energy barrier between these two regions is almost identical. For comparison, we also plot 
the tree-level effective potential with the running constants calculated at the one-loop level 
in figure 13. We see that the tree-level potential barrier is lower by roughly two orders of 
magnitude with respect to the one-loop case. Moreover, the instability region for the tree- 
level potential starts around h = 1.5 • 10 10 GeV and approximately coincides with the point 
at which A h becomes negative. On the other hand, for the one-loop potential this region 
is shifted towards the larger field value, namely h ~ 4.5 • 10 10 GeV. A similar conclusion 
concerning the influence of the higher loop corrections on the stability of the Higgs effective 
potential were obtained for the case of the Standard Model Higgs in flat spacetime [4]. 

5 Summary 

In this paper we have investigated the problem of the influence of the gravitational field 
on the stability of the Higgs one-loop effective potential. We focused on the effect of the 
classical curved background as opposed to the usual flat (Minkowski) background plus 
gravitons corrections. To this end, we used a local version of the heat kernel method, as 
introduced by DeWitt and Schwinger, which allows to investigate the case of large but slowly 
varying curvature of spacetime. To represent our quantum matter sector we used gauge-less 
top-Higgs sector (we chose the unitary gauge for the Higgs field and specialized to its radial 
mode). We also considered the presence of the second heavy real scalar coupled to the 
Higgs field via the quartic term. This scalar, when not possessing the vacuum expectation 
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Figure 12: The one-loop effective potential along the trajectory connecting the electroweak 
minimum and the region of the instability at high fields values. The running energy scale 
was set as p, = ^=/i. The spacetime curvature was given by the energy density p = (t^ 4 + ^ 4 , 
where a = 50 and v = 10 9 GeV. 



Figure 13: The one-loop (1A 1 ') and the tree level (1A°1) effective potentials along the 
trajectory connecting the electroweak minimum and the region of the instability at high 
fields values. The running energy scale was set as p = -^/t. The spacetime curvature was 
given by the energy density p = <rz^ 4 + /r 4 , where a = 50 and v = 10 9 GeV. The insert 
shows the behavior of potentials around the maximum of the tree-level potential. 
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value, may be dark matter candidate or when it possesses the vev it may be considered 
as the mediator to the dark matter sector. We focused on the latter case. Moreover, we 
considered both fields to be non-minimally coupled to gravity. 

Applying the heat kernel method, we obtained the divergent and finite (up to terms of 
the second order in curvatures) parts of the one-loop effective action. From the divergent 
part we got the beta functions for the theory at hand. We have found that, in agreement 
with the general results, the beta functions for various scalar quartic couplings, top Yukawa 
coupling and gamma functions for the scalars masses and field strength renormalization 
factors are the same as in the flat spacetime case. This is due to the fact that we considered 
purely classical gravitational background (without gravitons). We have also found beta 
functions for the non-minimal coupling constants [£,h/x) of the scalar fields in the model 
(3.19), (3.20). After investigating the running of these coupling constants we conclude that 
if we assume that they are initially zero (£,h/X ( m t) = 0, where mt is top mass) they run 
towards negative values at the high energy scale (figure 4b). On the other hand, if we 
postulate that they are initially above conformal value (^ = g) they run towards larger 
positive values in the high energy region (figure 4c). 

We have also given the explicit form of the one-loop effective action containing terms up 
to second order in curvatures. Namely, our action contains terms linear in the Ricci scalar 
(. R ), quadratic in the Ricci scalar and the Ricci tensor (R^) and linear in the Kretschmann 
scalar (/C = Rfu/a^R 111 ' 01 ^) (2.32). 

After confirming that, like in the flat spacetime case, our model possesses an insta¬ 
bility region for the large Higgs field value (figure 5b), we turned to the investigation of 
the influence of the gravity induced terms on the shape and the stability of the effective 
potential. Firstly, we considered the radiation dominated era and found that the one-loop 
induced terms (the tree-level ones are absent as the consequence of Friedman equations and 
the equation of state, namely in this era we have R = 0) give small positive contribution 
to the effective potential at the electroweak minimum (figure 6). The magnitude of this 
contribution is dependent on the total energy density. Figures 7 and 8 represent the same 
kind of effect but also for the de Sitter and matter dominated eras. The main difference 
between these two figures is the fact that for the first one we have £h/x( m t) = 0, while for 
the second one £h/x( m t) = g- In the absence of the tree-level terms {^h/X = 0 case) the 
gravitational terms contribute negatively to the effective potential in the de Sitter and mat¬ 
ter dominated eras. Moreover, this effect appears at the one-loop level. On the other hand, 
when £h/x >0 the gravity induced contributions are positive also for the aforementioned 
eras and they are in fact orders of magnitude bigger than for the radiation dominated era 
even for small values of X (£h/x 0(1)). 

The last problem relevant for the small field region which we considered was the influ¬ 
ence of the gravity induced terms on the critical temperature needed for the destruction of 
the electroweak minimum. Focusing on the qualitative description of the problem we have 
found the formulae for the critical temperature for the de Sitter and radiation dominated 
phases of the Universe evolution. They are given by expressions (4.28),(4.29) and (4.33), 
respectively. The obtained relations indicate that there are two types of corrections. The 
first one is additive and is suppressed by negative powers of the Planck mass. The sec- 
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ond one is multiplicative and is inversely proportional to the scalar non-minimal coupling 
constant (£h)- This type of correction is important for the de Sitter era and may change 
the critical temperature even by an order of magnitude (for large £) in comparison to the 
flat spacetime one. On the other hand, for the radiation dominated era we have only an 
additive negative contribution that is suppressed by Mp 4 . 

Since we used the truncated series representation of the heat kernel, a comment about 
the validity of presented results is in order. In fact, all the results summarized so far are 
obtained in the region where R < m 2 H , or 7Z 2 < for the radiation dominated era, 

where rn? H is the physical Higgs mass squared (m#_ ~ 125 GeV) and TZ 2 represents terms 
that are quadratic in Riemann and Ricci tensors. In this region our approximation is a very 
good one. 

We also pursued the question of how big energy density should be in order to in¬ 
duce a qualitative change in the one-loop effective potential for the scalar fields. To this 
end, we investigated regions of small (around electroweak minimum) and large (around 
instability scale) fields. In the small fields region we found that the gravity induced term 
contributes positively to the effective scalar mass parameters ( m(h)heff and m[h)x e ff) hi 
the Lagrangian if we are in the radiation dominated era or if we have a positive value of 
the non-minimal coupling constants in de Sitter and matter dominated eras. We defined 
the effective mass parameter in a manner similar to the definition of the effective quartic 
coupling in large field region, namely m(h ) 2 ^ = — Our calculations revealed that 

for the energy scale of the order v ~ 10 11 GeV, with the standard assumption that p = zz , 
this contribution is large enough to change the sign of rn {h)\/x e ff i which leads to the dis¬ 
appearance of the electroweak minimum. Since this energy scale lies slightly above the one 
allowed by our approximation (v ~ 10 9 GeV), we treat this result rather as an indication 
that gravity induced effect should be investigated more carefully even for the energy scales 
well below the Planck one than the statement of the actual effect. 

As far as the large held region is concerned, we investigated the influence of gravitational 
terms on the effective scalar quartic self-coupling of the Higgs held (dehned as A (h)h e ff = 
AV< h i ^ ' We P resen t e d results for the radiation dominated and de Sitter eras in figure 10 
and figure 11 . We found that for the sufficiently high energy density we get an improvement 
of the stability for the radiation dominated era and also for the de Sitter era for the positive 
non-minimal coupling constants. This means that gravity induced terms contribute positive 
factors to A (h)h e ff- On the other hand, if ^ is negative at large energy then the stability is 
worsened. We calculated the order of magnitude of the energy density for this effect to take 
place and we found that it is equivalent to the energy scale v r*j 10 13 410 14 GeV, while the 
Higgs held is of the order h ~ 10 10 GeV. This means that most energy is not stored in the 
Higgs held. Again, this is the above region of validity of our approximation v r-*j 10 10 GeV 
and should rather be treated as an indication of the possible effects. Nevertheless, we found 
it interesting that gravity may induce non-negligible effects at energy densities much below 
the Planck density, in the considered case we have p ~ 10 -21 4 lCV 20 /Jp, where pp is the 
Planck energy density. 

As the hnal remark we point out that it would be very interesting and important for the 
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problem of the stability of the Standard Model to go beyond limits of our approximation. 
Unfortunately, this requires another representation or a resummation technique of the heat 
kernel that could be applied to the case of large and slowly varying background fields, which 
at the present time we are unaware of. 
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